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For a class of analytical parametrizations of the freeze- 
out state of relativistic heavy ion collisions, we perform a 
simultaneous analysis of the single-particle mx-spectra and 
two-particle Bose-Einstein correlations measured in central 
Pb+Pb collisions at the CERN SPS. The analysis includes a 
full model parameter scan with \ 2 confidence levels. A com- 
parison of different transverse density profiles for the particle 
emission region allows for a quantitative discussion of possible 
model dependencies of the results. Our fit results suggest a 
low thermal freeze-out temperature T ~ 95 ± 15 MeV and a 
large average transverse flow velocity v± ~ 0.55 ±0.07. More- 
over, the fit favours a box-shaped transverse density profile 
over a Gaussian one. We discuss the origins and the conse- 
quences of these results in detail. In order to reproduce the 
measured pion multiplicity our model requires a positive pion 
chemical potential. A study of the pion phase-space density 
indicates ^ w 60 MeV for T = 100 MeV. 



I. INTRODUCTION 

The space-time analysis of hadronic one- and two- 
particle spectra measured in relativistic heavy ion col- 
lisions has attracted much attention in recent years. Via 
the reconstruction of the hadronic phase-space distribu- 
tion at freeze-out, this method can provide detailed geo- 
metrical and dynamical information about the last stage 
of the collision. The ultimate goal of the experimental 
relativistic heavy ion program is to produce and test the 
dense early stage of the collision in which quarks and 
gluons are the relevant degrees of freedom and a quark- 
gluon-plasma is expected (for an up-to-date overview see 
[ ^ |H). However, only the hadronized remnants of this 
state are experimentally accessible. Characterizing their 
spatial and dynamical distribution, a particle interferom- 
etry based space-time analysis can provide estimates of 
the phase-space density attained in the collision and it 
can establish an experimentally determined endpoint for 
microscopic simulations of the complicated multiparticle 
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dynamics. This makes it a valuable tool in the search for 
the quark-gluon-plasma. 

One of the main motivations for space-time analyses in 
recent years is the discovery that HBT (Hanbury Brown 
and Twiss) particle interferometry allows to disentangle 
between random and directed dynamical components in 
the collision [ ||, [|, ^, ||, [?], || . Specifically, in the context 
of hydrodynamic parametrizations which provide a very 
convenient characterization of the freeze-out region, the 
one-particle spectra are determined by an effective blue- 
shifted temperature from which temperature and flow ef- 
fects cannot be separated unambiguously [ ||]. On the 
other hand, the M^-dependence of the transverse corre- 
lation radii increases with transverse flow but decreases 
with stronger thermal smearing [ [l^] . By combining both 
observables the ambiguity between temperature and flow 
can be removed. 

Knowledge of the magnitude of the transverse flow is 
crucial for a dynamical picture of the transverse expan- 
sion of the collision system and for a dynamical back- 
extrapolation into the hot and dense early stage of the 
collision. With this motivation, there have been several 
recent discussions about disentangling temperature and 
flow effects. Except for the recent analysis in most of 
these discussions a re p ublished in conference proceedings 
and reviews [|?|, ||, Q |T3 , 14 1, and they are mainly based 
on preliminary data. None of them provides a full model 
parameter scan with \ 2 confidence levels for the two- 
particle spectra, and none of them makes quantitative 
statements about the model-dependence of the conclu- 
sions reached. However, both these points are very im- 
portant, since one may wonder, e.g., to what extent the 
main conclusions, a relatively low temperature and high 
transverse flow velocity, are subject to model-dependent 
details of the parametrization. The present work ad- 
dresses this gap in the literature with an extensive model 
study. 

Three collaborations measured and published data on 
Bose-Einstein correlations for 158 AGeV/c Pb+Pb col- 
lisions at the CERN SPS: WA98 [ fL§, NA44 I 01 and 
NA49 [ [llj]. Due to their small acceptance, NA44 can 
only determine particle correlations in two transverse mo- 
mentum bins which moreover correspond to slightly dif- 
ferent rapidities. Such data are not well suited for a de- 
tailed analysis where the Mj_-dependence of the correla- 
tion radii plays a crucial role. The situation is somewhat 
better for the recently presented WA98 data [ [l5| but 
these were not yet available at the time of the present 
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analysis. Most suited for our analysis tool is a large 
acceptance experiment like NA49 which covers almost 
the whole forward rapidity region. Their published data 
are, however, given only in the YKP parametrization [ 
jlTf . Although in many situations its parameters have 
the most straightforward physical interpretation [ p"7| , 
this parametrization can be ill-defined in some kinematic 
regions [ [l^, |l9) , and a cross-check with Cartesian (Pratt- 
Bertsch) correlation radii is hence necessary. The NA49 
experiment has the unique capability of cross-checking 
experimental results from different detector components 
with overlapping acceptance. At the moment, the re- 
maining differences between the different detectors are 
still larger than the published error bars [ |2(J. We ad- 
dress these subtleties in Section || where we discuss in 
detail the data used in our analysis. 

The class of fireball models used for our analysis was 
described in detail elsewhere (for recent reviews see e.g. 
[ [l^, [b|]). To get an idea of the degree of model depen- 
dence of our final conclusions we here investigate two dif- 
ferent transverse density distributions, a box-shaped and 
a Gaussian profile. The relevant model parameters are 
shortly introduced in Section HI. Up to now our analy- 



sis is restricted to the transverse momentum dependence 
in the particular rapidity bin (3.9 < Y < 4.4) for which 
a complete set of HBT parameters was published in [ 
ITU . A realistic description of the rapidity dependence 
of the transverse one^particle spectra would require a re- 
fined model f ||, |2l], |22[|. The results of the fit, as well as 



related technical details, are given in Section IV 



As an application of these results we discuss in Sec- 
tion |\]thc average pion phase-space density at freeze-out. 
Bcrtsch [ pointed out that knowledge of the single- 
particle spectrum and Bose-Einstein correlations allows 
for a model-independent extraction of this quantity. The 
reason is that absolutely normalized single-particle spec- 
tra carry information about the particle density in mo- 
mentum space while the width of the distribution in con- 
figuration space can be extracted from correlation stud- 
ies. We extract this model-independent quantity from 
the data and compare it with the prediction of our model; 
in this way we determine the value of the pion chem- 
ical potential needed to reproduce the measured pion 
yields. This is an interesting quantity since a large value 
of fi n would indicate an "overpopulation" of phase-space 
and could indicate the onset of multiparticlc effects (e.g. 
stimulated emission or a pion laser). 

Section VI contains our conclusions. Technical details 



are deferred to three Appendices. 



II. THE DATA 

Data on the Bose-Einstein correlation radii as func- 
tions of K± from the 5% most central Pb+Pb collisions 
at 158 A GeV/c were published by the NA49 collabo- 
ration [ |ll| for the rapidity window Y E (3.9, 4.4) (i.e. 



1 < Y CM < 1.5). Our analysis will focus on these data. 

The single-particle p^-distributions of negatively char- 
ged hadrons (h~) are taken from [Q. At the time of our 
analysis, spectra of identified pions were not yet available 
for the above rapidity window. The differences between 
h~ and identified pion spectra (from negative kaons and 
a few antiprotons) are, however, small since most nega- 
tively charged hadrons are pions. This is good since these 
differences will have to be modelled and thus introduce 
(small) systematic uncertainties. The rapidity binning 
of the h~ data in [ Q is slightly different from that of 
the correlation analysis of [ [TTJ . For h~ we use the bin 
4.15 < y < 4.65 which has the largest overlap with that 
of the correlation data. This is an excellent approxima- 
tion since the p \ -spectra vary only weakly with rapidity 
in this region [ |24j . 

For wide bins the question arises where exactly the 
data points should be placed when comparing them to 
model predictions. Usually one puts them in the middle 
of the bins. This may, however, lead to systematic errors [ 
|25| if the measured distribution is not flat. Our procedure 
"where to stick the data point" is state of the art [ p5| : 
if the measured distribution is well parameterized by a 
function g(x), then the appropriate position of the data 
point corresponding to a bin of width Ax between x\ and 
x-i is obtained as 
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Ax 



g(x) dx 



(2.1) 



where g~ x is the inverse function of g. The position of 
the data points of the single-particle spectrum in rapidity 
is calcul ated from this equation using for g{x) the distri- 
bution ( [All ) witri Ay = 1.4 [ ||, |27|; to get the appropri- 
ate p osit ions in transverse momentum, the parametriza- 
tion ([AJ) with T inv = 185 MeV (inferred from the mea- 
sured average transverse momentum and the fact that 
the observed pj_-spectrum is very well parameterized by 
an exponential with inverse slope Tj nv ) is used for g(x). 
Note that the width of the p± bins in the used data was 
100 MeV/c, with the first bin starting at 50 MeV/c. For 
these data we have explicitly checked that taking the data 
points in the centres of the bins or according to the above 
described procedure does not lead to an observable dif- 
ference in the fit results. 

The HBT two-particle correlation radii published in [ 
|IT| for the rapidity window 1 < Y CM < 1.5 have sev- 
eral shortcomings which can only be resolved in future 
experimental analyses: First, the correlation radii are 
only given for the YKP parametrization. Their statis- 
tical uncertainties are systematically larger than for the 
Cartesian parametrization [ p8j, ^9| , and the data shown 
do not allow for the important cross-check [ 
with the Cartesian parameterization. This is 
an important issue since the YKP fit was recently found 
to be problematic [ |TJ, ||. Second, Ref. [ |ll) contains 
only data taken in the Main Time Projection Chamber 
(MTPC), in spite of the observed systematic differences 
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between radii extracted from the MTPC [ p8[ and the 
VTPC data [ || (of the order of 0.5 fm [ |g^The error 
bars in [ do not include this systematic uncertainty 
which is only roughly estimated to be about 15%. 

To address these difficulties we based our analysis on 
the three accessible, but so far unpublished complete data 
sets from the NA49 experiment: h + h + and h~h~ corre- 
lations from the VTPC [|29| and hrhr correlations from 
the MTPC [ The reported errors are the output of 
the MINUIT fitting routine which is known to underes- 
timate errors [ ||(J . For our analysis we considered these 
three data sets as independent measurements and took 
their average, with correspondingly increased error bars 
which now also include the systematic deviations between 
these measurement. Since both in theory and experiment 
the squared correlation radii and the YK velocity are the 
directly determined fit parameters [ [l7], |3l| , we average 
over and fit to those. 

The obtained fit parameters for the Cartes ian and 
YKP parametrizations are summarized in Tables 



C2, respectively. They are also displayed in Figs 




The minor difference between the rapidity bins of the 
VTPC analysis [ |§ (3.9 < Y < 4.4) and the MTPC 
analysis [ g8| (4 < Y < 4.5) is irrelevant and has been 
neglected. The binning in transverse momentum differs 
slightly, too. Only four of the total number of five K± 
bins could be taken into account in the averaging proce- 
dure, and for the fourth bin a slig ht size difference be- 
tween [ and (see Table C2) had to be neglected. 
To position the data points in the bins we assumed here 
that their distribution g(x) inside the bin can be lin- 
earized. Then their position in the bin [xi,X2] can be 
calculated from [ E5| 



J^xg(x)dx 



dx 



(2.2) 



Inserting for t he f unction g(x) the di strib ution Q%(Y) 
from equation (A3) with Ay = 1.4 into (2.2), the position 
of the data points in the pair rapidity Y is obtained: 
V bin = 1.22. The values of AT 1 ?' 11 for the different K±- 

CM J- - ± - 

bins are calculated using for g{x) the function p^-(Aj_) 
from equation (A9), again with T; nv = 18 5 M eV. The 
resulting values are displayed in Tables C 2 and |C 2| . We 
observed [|l!| that slightly different (worse) model fits are 
obtained if the data points are placed at the bin centres 



rather than at the positions calculated from (2.2) 



It remains to check the compatibility of the Cartesian 
and YKP correlation parameters in order to ascertain 
the validity of the YKP data. To this end the Carte- 
sian radii were calculated from the YKP parameters and 
vice versa via cross-check relations published in [ |l7], [l8| . 
The resulting calculated radius parameters are shown in 
Figs. |l] and |^ by gray dashed symbols. The errors were 
propagated to the calculated radii via 
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(2.3) 



where at is the calculated error of the Cartesian (YKP) 
correlation parameter R 2 and a'j the "measured" error 
of the YKP (Cartesian) parameter R' 2 . (Here the notion 
'YKP correlation parameter' includes radius parameters 
as well as the YK velocity v YK .) Note that the non- 
diagonal terms of the error matrix are not known and 
had to be neglected. In most cases the calculated pa- 
rameters show larger errors than the directly measured 
ones. This might improve if non-diagonal error matrix 
elements could be taken into account in Eq. (2.3). 

Figs. [I] and || show that the YKP and Cartesian data 



are "on average" consistent [ 32 . In view of the known 
fragility of YKP fits this is an important and non-trivial 
result. However, it was recently argued [|l^, |l^| that the 
YKP parametrization is quite subtle and can become ill- 
defined for certain (even realistic) sources. The only way 
to test this possibility in experiment is to calculate the 
YKP parameters from the Cartesian ones and check that 
the resulting values are real. Since experimental Carte- 
sian radius parameters have a finite measurement error, it 
is not sufficient to perform this consistency check only for 
their average values; rather, all parameter values within 
the error interval should be checked. When doing so we 
found problems with the definition of the YKP parame- 
ters inside the error intervals for the third and fourth K± 
bins. These may be related to our further observation [ 
[l9| that no good model fit was possible starting from the 
YKP parameters, and that the fits tended to drift into 
strange parameter regions. This emphasizes the future 
need for an explicit cross-check between Cartesian and 
YKP parameter fits directly on the experimental level, 
including error propagation with the complete error ma- 
trix. Due to our problems with the measured YKP pa- 
rameters, the model fits presented here are based on an 
analysis of the Cartesian correlation radii. 



III. THE MODEL 

A. The emission function 

Our model analysis is based on the widely used emis- 
sion function [ § [§ H H, g H 

Six, K) d±x = M±C °Sf-^ exp (- K ^ (X) + | 

{V - Vo) 2 



2{Ar 1 ) 2 



T 

drj r dr dip 



(2tt) 3 
xG(r) exp 

*ot™H-wJ' (3 ' 1) 

Here, the pair momentum K is parameterized in terms 
of the transverse mass M± = -jK'j_ + m 2 , longitudinal 
rapidity Y, and the transverse momentum K± 
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K = (Mi coshF, K x , 0, Mi sinhF) . 



(3.2) 



We use here the usual coordinate system with x-axis di- 
rected parallel to K j_. In configuration space, we use the 
polar coordinates r and tp for the transverse directions, 
a longitudinal proper time r = \Ji 2 — z 2 and space-time 
rapidity 77 = | In . 

The 77-profilc in the model is chosen to be a Gaus- 
sian peaked at 770 with width Ar/. In this way 770 de- 
termines the centre of mass (CMS) rapidity. Since we 
focus on a single rapidity bin and work in the LCMS, 
we set 770 = —1.22 (see previous Section) and A?7 = 1.3. 
The latter value is obtained from a comparison with the 
single-particle rapidity distribution. 

Freeze-out proper times are distributed according to 
a Gaussian of width At centred at tq. This allows for 
freeze-out of particles during an extended time period 
given by At. 

The transv erse geometry is specified by the distribu- 
tion G{r) 



As only azimuthally symmetric fire- 
balls are considered, this function does not depend on p. 
In our analysis two particular profiles will be assumed: a 
Gaussian one 



G(r) = exp 
and a box-shaped one 



2R% 



G(r) = 9(R B - r) . 



(3.3) 



(3.4) 



The fireball at freeze-out is assumed to be in local ther- 
mal equilibrium at temp erat ure T. This is modelled by 
the Boltzmann factor in (3.1). The chemical potential /1 
in its argument has no effect on the correlation radii and 
only affects the normalization of the single particle spec- 
trum. The argument K-u(x) allows for longitudinal and 
transverse expansion of the fireball through the collective 
four- velocity field u(x): 



u(x) — (cosh 774 cosh 77, cost/? sinh?7 t , simp sinh?7 f 
cosh?7 t sinh?7) . 



(3.5) 



The longitudinal expansion rapidity tji (x) has here been 
identified with the space-time rapidity 77. This leads to 
a longitudinal expansion velocity vi ong (z,t) — tanh.77 = 
z/t, corresponding to Bjorken-type boost-invariant ex- 
pansion [ B6| . The transverse expansion is parameterized 
by the transverse flow rapidity r) t (x). It is assumed to in- 
crease linearly with the distance from the collision axis, 



(3.6) 



The scaling factor 77/ specifies the value of the transverse 
flow rapidity at the transverse rms radius, given by 



r rms = V2R G 
for the Gaussian transverse distribution and by 



(3.7) 



Rb 
75 



(3. 



for the box-shaped one. Fr om these relat ions we expect 
the two fit parameters in ( |3.3| ) and (3.4) to satisfy ap- 
proximately Rb ~ 2 Ra- 
in the literature the transverse flow is often quoted in 
terms of the average transverse expansion velocity v± 



/ °° r dr tanh r\ t (r) G(r) 

/-no — , , : 



Jo r 



drG(r) 



(3.9) 



Its value is roughly given by the radial velocity at 
the transverse rms radius, v±(r lms ) ~ tanh^j; look- 
ing more closely, it is slightly smaller: for 7/^=0.6 one 
finds i>_L(r rms )«0.54 and vj_=0.50 (0.46) for the box-like 
(Gaussian) transverse density profile. 

Our calculation of the single-particle spectra from 
( |3.1[ ) includes contributions from resonance decays as de- 
scribed in [ ^|, H?| . This is crucial for a correct description 
of the yield and shape of the spectra. We include all de- 
cays with branching ratios above 1% from mesons with 
masses up to 1020 MeV/c 2 and from baryons with masses 
up to 1400 MeV/c 2 . For decay chains, the product of the 
branching ratios is required to be larger than 1%. For 
baryons the ch emi cal potential [1b {T) was parameterized 
according to (|Bl|). Strange particles acquire a chemi- 
cal potential ns(T) determined by the condition (B2) of 
strangeness neutrality. Finally, contributions from nega- 
tive kaons and antiprotons were included in the calcula- 
tion of the h~ spectrum. 

The role of resonance decay contributions to the HBT 
radius parameters is known to be much less important [ 
|37| . Our calculation of the correlation radii will thus only 
include direct pions. This is an essential technical simpli- 
fication: the additional integrals from the resonance de- 
cay phase space would have increased the numerical task 
from calculating a 2-dimensional integral to calculating a 
5-dimensional (for two-particle decays) or 6-dimensional 
one (for three-particle decays) . This very time consuming 
calculation was performed in model studies [ ^ where 
the correlator was evaluated for only a few characteris- 
tic sets of model parameters. However, in a simultaneous 
multi-parameter fit, in which the fit routine calls the two- 
particle correlator approximately 50-100 times for thou- 
sands of different model parameter combinations, it can- 
not be done. 



We finally remark that the model (3.1) is customarily 
used with a Boltzmann distribution rather than a Bosc- 
Einstein one, since this greatly simplifies the analytical 
and numerical analysis. The resulting differences are usu- 
ally small. If the pions develop a positive chemical poten- 
tial, however, this visibly affects the shape of the single 
particle spectra; for this reason we use, in fact, a Bose- 
Einstein distribution to calculate the pion multiplicities 
and the single-particle spectra of direct pions (see Ap- 
pendix [b] and |c]) . For the heavier resonances the Boltz- 
mann approximation is excellent. For technical reasons 
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we also use the Boltzmann approximation in the calcula- 
tion of the correlation radii, in line with the much larger 
systematic error of the corresponding available data. 



B. Basic relations 

We shortly recall how the emission function S(x, K) 
is related with the observables to be calculated. More 
details can be found e.g. in [ fL2|, 

The single-particle spectrum is obtained from 



dN 



d x S(x,p) . 



(3.10) 



For the two-particle correlations, we used the Cartesian 
parametrization: 



C(q,K) - 1 = cxp [-q 2 Rl(K) - q 2 s R 2 s (K) 

-qfRf(K)-2q oqi R 2 ol (K)] 



(3.11) 



where the qi are the components of the momentum differ- 
ence in the out-side-long coordinate system and K stands 
for the average pair momentum. The correlation radii are 
obtained from the emission function via 



if) 

((£■ 

((2- 



■ PJ) 2 ) , 
Pit) 2 ) , 
-Ait) (5 



Pit)) 



Here 



and 



(f(x))(K) = 



/ d 4 xf(x) S(x,K) 
jd A x S(x,K) 



(3.12a) 
(3.12b) 
(3.12c) 
(3.12d) 



(3.13) 



(3.14) 



In our calculations all integrations were performed nu- 
merically. This distinguishes our work from Ref. [ Q 
where analytic approximations were used instead. For 
sources with strong transverse flow (as will be the case 
here) the analytical approximations for the correlation 
radii may become problematic [ |3^, |34|, Also, it is 
not easy to capture the intricate effects of resonance de- 
cay kinematics on the single-particle spectra with simple 
analytical expressions as those used in [ . 



IV. THE FIT 

A. Single-particle px-spectrum 

The first step is to fit the h~ single-particle p±- 
spectrum with the model (3.1). Resonance decay con- 
tributions are treated as described in Sec. Ill A. The 



geometrical model parameters A77, Rb (Rg)i tq, At are 
known to only affect the absolute normalization of the 
single-particle spectrum, but not its shape [ ||, [l^] . For 
fixed transverse density and flow profiles, the spectral 
shape is completely determined by the temperature T 
and the transverse flow strength rjf. 

The fit to the spectrum was performed with the CERN 
package MINUIT [ 39 , using for the calculation of ( 3.10| ) 
(see Appendix [b]) a routine described in [ [37). The 
data were compared with the y-integrated m^-spectrum 
rather than with the spectrum at a fixed value of y. Since 
the y-integration can be done analytically and removes 
the dependence on the longitudinal kinematic limits for 
resonance decays, the former requires much less computa- 
tion time. The differenc es a re very small since the m±- 
spectra obtained from (3.1) are almost y-independent, 
except for very forward/backward rapidities. 

The resulting \ 2 contour plots in T and rjf are shown 
in Figs. |^ (box-shaped transverse geometric profile) and 
^ (Gaussian profile). For comparison the average trans- 
verse expansion velocity v± is given on the right ordi- 
nates of these Figures. For their calculation the routine 
was driven through the whole T-rjf domain covered by 
the Figures, and in each point the normalization of the 
spectrum was fitted as a third parameter. One finds a 
clear x 2 "valley" pointing from the upper left to the lower 
right corner, reflecting the anticorrelation between T and 
rjf in the m^-slope. The indicated confidence levels re- 
sult from a ^-distribution with 16 degrees of freedom 
(19 data points minus 3 fit parameters) 
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The two panels of Fig. ||| show the effects of including 
the K~ and p contributions in the negative hadron spec- 
tra - the upper panel includes only negative pions (in- 
cluding all resonance decays). At low temperature the fit 
results are nearly identical because the K~ and p con- 
tributions are strongly suppressed by their masses and 
by the large baryon chemical potential (at T=80 MeV 
we have jiB w 430 MeV), but at higher temperatures 
their inclusion reduces appreciably the amount of trans- 
verse flow needed to fit the measured slope. Our re- 
sults differ in two points from those published in [ 
(i) We include nonzero baryon and strangeness chemical 
potentials. This increases the decay contributions from 
baryon resonances, making the spectra steeper and thus 
requiring more transverse flow to reproduce the measured 
slope. Compared to the case of vanishing chemical po- 
tentials this shifts the "% 2 valley" in Fig. |3|a upwards by 
« 0.05 c (somewhat less at low T and more for higher 
T). (ii) We include contributions from K~ and p. This 
flattens the //"-spectrum further because these heavier 
particles are more strongly affected by transverse flow. 
To fit a given spectral slope, lower values of rjf are thus 
needed at given T. The two effects (i) and (ii) are seen 
to partially cancel each other. 

A comparison of Figs. I and | gives a feeling for the 
model-dependence of the fit resulting from two different 
choices for the transverse density profile. Since at fixed 
rjf they lead to different average transverse flow veloci- 
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ties v± (the quantity which determines the blueshift of 
the spectral slope), one should use the labelling on the 
r.h.s. of these Figures for comparison. One sees that for 
a Gaussian density profile slightly smaller values for v± 
are needed to account for the measured slope. The rea- 
son is that superimposing a Gaussian distribution G(r) 
with a monotonically increasing transverse flow profile, 
a significant part of the high-pj^ part of the spectrum is 
obtained from contributions at large transverse distances 
r > Rq in the Gaussian tails. A box profile does not 
allow for emission from distances r > Rb and thus re- 
quires slightly larger values for v± to account for the same 
spectra. 

Some access to the shape of the transverse density 
profi le is possible through two-particle correlations (see 
Sec. tVB ), but limited statistics and other uncertainties 
leave some room for interpretation. The differences be- 
tween Figs. | and | should thus be mainly taken as an 
estimate for the systematic model uncertainties in the fi- 
nally extracted values for T and v± . That they are small 
for the single-particle fits is certainly gratifying. 



B. Two-particle correlations 



The correlation radii summarized in Table C 2 were fit- 
ted with the CERN package MINUIT [ |§, using for the 
model calculatio n a r outine which computes the correla- 
tion radii from ( 3.12 ). We scanned the whole T-rjf do- 
main as before, performing in each point a 3-parameter fit 
to find the best values of Rb (Rg), r o, and At. A77 = 1.3 



was always kept fixed, see Sec. [II. 



The results of these fits, superimposed on the fit of the 
single-particle spectra, are shown in Figs. || (box- profile) 
and ^ (Gaussian profile). The contours correspond to a 
X 2 -distribution of 11 degrees of freedom (16 data points 
minus 5 model param eters) . These results will be fur- 
ther discussed in Sec. IV D. Here we only observe that 



the HBT radii indeed allow to disentangle the ambiguity 
between temperature and transverse flow (although the 
uncertainties are still significant). The box-shaped trans- 
verse geometric profile seems to be favoured by the fit, 
although limited statistics does not allow to rule out a 
Gaussian shape. Independent of the choice of the trans- 
verse density profile we can, however, safely conclude that 
the data require strong transverse flow with v± > 0.3. 
The best fits favour low freeze-out temperatures between 
80 and 110 MeV and large average transverse expansion 
velocities between 0.47c to 0.62c (for the box model). 

On the other hand, the "x 2 valley" defined by the cor- 
relation radii is clearly seen to deviate from the simple 
dependence on rjf and T which one obtains in analytical 
approximation for Gaussian transverse density profiles [ 

H0: 



R 2 s (Mx) 



R 2 G 



2T' 



Mi 



(4.1) 



Analyses employing this relation (e.g. [ |l4], pT| ) should 
therefore be taken with some caution. 



C. Total yield 

Our fit determines for each combination (T, rjf) of tem- 
perature and transverse flow the remaining model para- 
meters. This allows for the calculation of the total pion 
multiplicity by integrating the model emission function 
over the whole phase-space. The corresponding formu- 
lae are derived in Appendix |c[ We recall that for this 
calculation we use a Bose-Einstein distribution for the 
direct pions and that resonance decay contributions are 
included. 

The resulting 7r~ multiplicities (including pions from 
resonance decays) are shown in Fig. [7] for both box-like 
and Gaussian transverse density profiles. At very low 
temperatures (up to 100 MeV) the total pion multiplicity 
is dominated by directly produced pions. With increasing 
temperature the fraction of pions from resonance decays 
grows rapidly. For the model with a box-shaped density 
profile the multiplicity grows faster with temperature; for 
rjf = this model produces (in Boltzmann approxima- 
tion) twice as many pions as a Gaussian model whith 
the same transverse rms radius (i.e. with Rq = Rb/2). 
This is due to the larger covariant volume occupied by 
the box-shaped model. 

The experimentally measured negative hadron multi- 
plicity 715 ± 30 [ |2?J is indicated by thick contour 
lines. Clearly, these "bands" lie far outside the regions of 
(T,rjf) favoured by the fit to the spectrum and correla- 
tion radii (cf. Figs. ^| and [|). They tend to favour much 
higher temperatures. Again, the model with a box-like 
density profile is favored since it minimizes this discre- 
pancy. Nevertheless, a non-zero pion chemical potential 
must be introduced in both models for a correct repro- 
duction of the measured multiplicity. This question is 
studied in more detail in Sec. [v|. 



D. Results and discussion 

In order to obtain a more direct picture of the quality 
of the above fits, we selected three sets of fit parame- 
ters for each of the two models (Gaussian (gl, g2, g3) 
and box-shaped (bl, b2, b3)), indicated in Figs. |B| and[(j. 
The c orre sponding complete parameter sets are given in 
Table ID which also shows the predicted total ir~ multi- 
plicities (including resonance decays) in each case. 

The quality of the fits bl and gl to the measured single- 
particle m^-spectrum [ ^4| is seen in Fig. ||. Similarly, 
Figs. [| (box) and |Io| (Gauss) show the quality of the 
fits to the Cartesian correlation radii. The box-shaped 
density profile accommodates the Mj_-dependence of R Q 
and R s better than the Gaussian one. It reproduces the 
rapid initial increase of R a , its rather steep decrease at 
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larger M± , and the slope of R s all reasonably well. With 
the (T, rjf) values allowed by the single-particle spectrum, 
the Gaussian model has difficulties in reproducing the 
strong M^-dependence in R s and R . 

The different fit qualities achieved by these two models 
reflect their different behaviour in the presence of trans- 
verse flow. This is illustrated in Figure [If], where trans- 
verse cuts through the effective emission region for parti- 
cles with K±=500 MeV/c (in ^-direction) are shown. For 
the Gaussian source the effective emission region moves 
parallel to K j_ outward into the tail of the density pro- 
file; its "outward" homogeneity length is larger than its 
"sideward" one. For the box-shaped distribution, which 
forbids particle emission from r > Rb, the opposite is 
found: the effective emission region gets squeezed to- 
wards the edge of the box, and the "outward" homo- 
geneity length is now smaller than the "sideward" one. 
This causes a more rapid decrease of R with rising K± 
in this case, and it is this feature of the model which is 
preferred by the data. 

Recent studies of deuteron production via coalescence 
appear to support our conclusion. They also favour box- 
shaped transverse density distributions, although for a 
different reason [ ^2|, [f3|: the box-profile gives more 
weight to regions of large transverse flow velocities, and 
only in this way can one understand the observed flatten- 
ing of the deuteron m^-spectra compared to the proton 
ones [ [l2| ^3j. In contrast, what matters here for the 
correlation radii is the increasingly negative contribution 
from (x 2 — y 2 ) to R 2 — R 2 at large Mj_. This property 



was originally attributed in [ 44 to opacity effects in the 
source, i.e. to the suppression of particle emission from 
the interior of the source in favor of surface emission. In 
this sense, a radially expanding source with a box-like 
density profile looks at large K± like an "opaque source" 
(see Fig. 0) . 

This ambiguity illustrates that particle interferomet- 
ric measurements can make statements about the rms 
widths of the homogeneity regions, but that it is a model- 
dependent task to interpret how these homogeneity re- 
gions and their .ff -dependences are generated. That 
transverse density profiles with a sharper edge than the 
Gaussian one can naturally explain both the "opacity 
effects" in the correlation radii at large K± and the 
stronger flow experienced by the deuterons may be taken 
as an indication, but not as proof that they form indeed 
the preferred parametrization of the source at freeze-out. 

Before turning to a discussion of the remaining two ra- 
dius parameters, Ri and R 2 ^ let us stress another very 
important feature of the source: its strong transverse 
growth before freeze-out. Both the box- like and Gaus- 
sian transverse density profiles give at freeze-out trans- 
verse rms radii r rms = RgV% — Rb/V% of w 9 fm. This 
is about twice the transverse rms radius of the original 
overlapping Pb-nuclei of ss 4.5 fm. In view of the re- 
sults obtained in [ [37j for a Gaussian density profile it 
seems unlikely that the neglected resonance decay contri- 
butions can account for a significant fraction of this large 



difference, although we did not check this possibility ex- 
plicitly again for the box-like distribution. Dynamical 
consistency requires that such a strong geometric growth 
is accompanied by strong radial flow, as indeed seen in 
our analysis. 

The longitudinal radius Ri is fitted very well by all 
selected parameter sets. This supports the validity of the 
assumed Bjorkcn scenario of boost-invariant longitudinal 
expansion of the reaction zone at thermal freeze-out. The 
measurement of the time parameters of the model (tq 
and At) is, however, affected by rather large errors, and 
it is model-dependent. The model-dependence of At was 
extensively discussed in [ [l8| to which we refer for details. 
Here we concentrate on a discussion of the parameter To- 

For a system undergoing boost-invariant longitudinal 
expansion, the size of Ri is dominated by the longitudinal 
flow velocity gradient [ 35 whose inverse grows linearly 
with the longitudinal proper time r. Under these condi- 
tions it is suggestive to interpret tq (as extracted from 
Ri) as the total time from impact to freeze-out [ [u^. For 
the sets of fit parameters listed in Table III this interpre- 



tation runs, however, into trouble: taking the measured 
final average transverse flow velocity of v±_ 0.5 c and 
assuming constant acceleration one would expect [ ^] 
the mean (rms) radius of the matter to expand roughly 
according to r rms (r) = r rms (0) + §u_l(t)t. (This slightly 
exaggerates the point to be made since the higher pres- 
sure will lead initially to stronger acceleration.) This ex- 
pression should be evaluated at the average emission time 
which for the source (3.1) is given by (t) = tq + . 
For the two preferred parameter sets bl and gl in Ta- 
ble |U we get (t) rs 8.4fm/c. With r rms (0) w 4.5 fm one 
thus finds r rms ((r)) ss 6.6 fm; this falls clearly short of 
the measured rms radius of 8.5 — 9 fm. 

For realistic expansion scenarios the time necessary to 
expand to such large radii thus exceeds both to and (t). 
(Only if the mean transverse flow velocity of v±=0.5 c had 
been established directly after impact, the matter could 
have expanded transversely from 4.5 fm to about 9 fm 
during the time (r) ~ 8.4fm/c.) The important impli- 
cation is that, although the reaction zone appears to be 
expanding boost-invariantly at freeze-out, it cannot have 
expanded so rapidly throughout its history. The fireball 
must have undergone longitudinal acceleration. Again, 
this should not surprise anybody: the measured trans- 
verse flow must have been created by transverse pressure 
gradients, and since pressure is locally isotropic, it must 
have also pushed longitudinally. The longitudinal veloc- 
ity gradient measured by Ri is thus only a snap-shot at 
freeze-out, and it is very likely that some fraction of the 
longitudinal flow, like the transverse one, has developed 
gradually by work done by the pressure. As sketched in 
Fig. [f2] this automatically leads to a longer real fireball 
lifetime, r roa i > tq, between impact and freeze-out. 
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V. THE PHASE-SPACE DENSITY 

From the measured single-particle momentum spec- 
trum and two-particle correlation function one can in- 
fer the phase-space density, spatially averaged over the 
homogeneity volume, of the particles immediately after 
freeze-out [ ^3|. Hence, this measurement is not only 
sensitive to the absolute multiplicity of the particles, but 
gives also hints about the possible appearence of "over- 
population" in some parts of phase-space. Such an over- 
population might give rise to "pion laser" phenomena 
etc. [ 0. Since we have seen that the pionic phase- 
space seems to be populated more densely than expected 
in chemical equilibrium, an analysis of the phase-space 
density appears to be of interest. 

In this section we first elaborate on the formalism de- 
veloped in [ |2^] and then apply it to the data. A very 
qualitative study [ |47) of a wide set of different measure- 
ments did not indicate a large overpopulation of the pio- 
nic phase-space, but showed signatures of the transverse 
expansion in form of a "flattening" of the TOx-depcndence 
of the (position-) averaged phase-space density. Here we 
want to perform, for the specific set of data analyzed in 
this paper, a more quantitative study of both these phe- 
nomena. Unfortunately, the data quality does not yet 
allow for high precision investigations; we have to keep 
this in mind and will stay on a rather superficial level. 



A. Formalism 

Combining the expression for the two-particle correla- 
tion function 



C(p u p 2 ) = 1 + 



|J d 4 xe ic > x S(x,K)\ 
Jd*xS(x, Pl ) Jd*yS(y, P2 ) 



(5.1) 



(q = p\ — P 2, K = (pi + p 2 )/2) with that of the one- 
particle spectrum (3.1C), we obtain 



P 1 (p 1 )P 1 (p2)[C( Pl ,p 2 )-l} = 



J d 4 x e lq - x S(x, K) 



(5.2) 



Let us introduce the time integrated emission function 

/oo 
dtS(t,x + vt,K), (5.3) 
-00 

where v = K / K°. This quantity allows us to rewrite 
the integral of (5J2) over on-shell momenta q satisfying 
q-K = 0, 

J A 5(q ■ K) [Px{K + §) Px{K - §) (C(q, K) - 1)] 
n^j#xT?{x,K). (5.4) 



Here, the on-shell approximation (which is valid for 
q 2 <C 4Ejj) has allowed the replacement K° — > Ek = 
\Jvn? + K 2 . We also find 



P 1 (p) = J d 3 xX(x,p) 



(5.5) 



T,(x,p) is not the phase-space density. In the following 
we establish how it is connected to the latter. The phase- 
space density f(t,x,p) is obtained by summing over all 
particles of a given momentum emitted up to time t along 
a given trajectory: 



f(t,x,p) 



(2tt) 



d£ S(t',x + v(t' -t),p) . (5.6) 



The factor in front of the integral assures the correct 
normalization of / to the number of particles for t > tf, 
where tf is the last instant of the freeze-out process. 
One easily can show [ |l9) that 

J d 3 xr(t>t f ,x,p) = J d 3 x^(x,p). (5.7) 



From (5.4) and (5.5) then follows 



d 4 qS(q-K) [P 2 (K)(C(q,K)-l)] 
w 7S3 J d 3 xf 2 (t>t f ,x,K), 
Px{K) = ^ J d 3 xf(t>t f ,x,K). 



(5.8) 
(5.9) 



In ( p.8|) we have performed the smoothness approxima- 
tionPi(if+§) w P x {K-f) w P\{K). Dividing these two 
equations (and changing the notation K—>p) we obtain 

J d 3 xf 2 {t > t f ,x,p) 



(f)(p) 



J d 3 xf(t > t f ,x,p) 
Pi(p) / d 4 q6(q-p)[C(q,p)-l] 



(5.10) 



(5.11) 



This allows to determine the phase-space density of free- 
streaming particles averaged over positions at constant 
global time, since all quantities on the r.h.s. can be mea- 
sured. Due to Liouville's theorem, the phase-space den- 
sity of free -streaming particles does not change, and 
hence ( [5.11 ) gives the phase-space density averaged along 
the freeze-out hypersurface. 

Indeed, for a hypersurface uj on which the freeze-out 
process is just completed and a global time coordinate is 
tf(x) one can show [ [l9| that 

E p I d 3 xf n {t>t f ,x,p) 

P»d 3 o-^x)f n (t f (x) 7 x 7 p) 7 (5.12) 
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where at is the hypersurface given by t = const. > tf and 
d 3 a^ is the infinitesimal normal vector to cry. The factor 
p-d 3 a is known from the formalism of Cooper and Frye 
[ Q and stands for the flux of t he pa rticles across 07. 
This relation allows us to rewrite (pTG) as 



_ S af P-d 3 <r(z)f 2 (tf(x),x,p) 
U)[P) f af P-dMx)f(tf(*),x,p) 



(5.13) 



This is the phase-space density averaged over the hy- 
persurface along which the freeze-out is just completed. 
Eq. ( J5.13 ) establishes what can be learnt about the phase- 
space density in a model-independent way. A back- 
extrapolation across the freeze-out boundary is only pos- 
sible if additional assumptions about the mechanism of 
particle production are made. 



If the correlation functi on is parametrized as in (3.11) 
the integration over q in ( |5.11 ) is simple and leads to 



(f)(.p) = -^r — 



d 3 N 



E p p± dp± dy < 



(5.14) 



However, in this form all pions count towards (5.14) irre- 
spective of their origin. We want to eliminate the pions 
from longlived resonances since they do not contribute 
to the phase-space density at freeze-out. (We keep pions 
from short-lived resonances because these decay essen- 
tially in the same spatial region where the direct pions 
are set free.) The standard procedure for separating off 
long-lived resonances is based on the observation that 
they lower the intercept Adh- of the correlator at q = 
[ @, 37 . Assuming that there is no coherent contribu- 



tion to pion emission which would lower the intercept 
even further [ [l9|, we obtain the phase-space dens ity of 
"direct" pions by multiplying the r.h.s. of ( |5.14 ) with 
V Adir [ ^]. We further simplify ( 5.14 ) by introduc- 
ing the following parametrization for the single-particle 
spectrum: 



d 3 N 



1 



p± dp±dy< 



1 dN 1 



^ dy T?{y) 



exp 



in 



The final formula used in the following data analysis 
then reads 



</>(p) = ^ 



2 E p Tl v (y) 



exp 
1 



P± 

-L mv 



dN 
dy 



(5.16) 



1. Data choice 

The rapidity distribution dN/ dy is estimated from re- 
sults given in [ |2(| |27j to be 139 ± 22. An estimate for 
the inverse slope is obtained from the measured (p±) [ |24| ] 
as Ti nv = 18 5 ± 10 MeV. The correlation radii are taken 
Since for the investigated rapidity bin 



from Table C 2 



no intercept parameter was given for the MTPC anal- 
ysis [ we estimated it by averaging only the VTPC 
data for negative and positive hadrons [ |2{| . As seen in 
Table IV they are fairly i^-independent; moreover they 



agree well with the intercept parameters extracted in 
M for 4<F<5. 



2. Phase-space densities from statistical distributions and 
realistic emission functions 

We first study the simple question to what extent 
the data support the simple assumption that the ob- 
served phase-space density follows a purely quantum- 
statistical distribution. Rough agreement of the data 
with a Bose-Einstein distribution was observed in [ [47| . 
Since our model reproduces both the measured one- and 
two-particle spectra, it allows us to refine these observa- 
tions. To this aim, we compare in a first step in Fig. [l3| 
the data with simple expectations from statistical distri- 
butions. 

Since in LCMS the energy coincides with the transverse 
mass E — , the phase-space occupancy following from 
a Boltzmann distribution is simply 



(/b)(p) = exp 



T 



(5.17) 



while the expectation based on the Bose-Einstein distri- 
bution reads 



(IbeKp) = 



1 



e t — 1 



(5.18) 



The temperatures for the corresponding curves in Fig. |l3j 
are taken from the parameter sets indicated in the up- 
per right corner. One sees that unless the temperature 
reaches 160 MeV (which is disfavoured by the fits) the 
slope of the calculated distributions is steeper than that 
of the data. 



The results of Sec. IV identify as the main source 
of this discrepancy the strong transverse flow at freeze- 
out. That transverse flow has a strong effect on (f)(p) 
can b e seen already from a minimal modification of 
Q5.17 ), (5-18), in which the statistical distributions are 
"boosted" to the rest frame of the point of maximum 
cmissivity for a given flow field, it max : 



B. Application 



(/b'Xp) = exp (- 

(f B E>)(p) = 



p ■ u n 



T 



( p-Umt 

exp \ t 



(5.19) 
(5.20) 
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We note that the shape of the expectations (/b')(p) 
and (fBE')(p) (e.g. "the bump" ) shown in Fig. |l5| depend 
on the determination of u max which is model-dependent. 
For instance, the bump vanishes if a Gaussian transverse 
density profile is taken, while the slope shows a qual- 
itativel y s imilar behaviour [ In a ny ca se, as seen 
in Fig. [13|, the statistical distributions Cl9|) and fl||(] ) 
overpredict the phase-space density significantly for large 
transverse momenta. A more refined model is needed to 
account for the data. 



by recalling that the distributions (5.17) and (5.18) were 
evaluated with energy distributions at the point of max- 
imum emissivity. They thus overestimate the data by 
construction. To understand the rather large chemical 
potentials required for our model, we recall that in the 
present study the value of this chemical potential is sub- 
ject to other model-dependent features: 



This was done, using the model described in Sec. Ill to 
calculate the HBT radius parameters and the o ne-pa rticle 
spectra which enters the numerator of Eq. ( |5.14 ). As 
mentioned above, the data are multiplied with y/ Adh- (p) 
and measure the phase-space occupancy of direct pi- 
ous only. Furthermore, without chemical potentials, our 
model leads to smaller multiplicities than measured. This 
allows us to determine the chemical potential needed in 
order to reproduce the observed phase-space density. If a 
non-zero chemical potential is considered, the difference 
between the Bose-Einstein distribution and its Boltz- 
mann approximation increases. Therefore, for the calcu- 
lation of the single-particle spectrum the Bose-Einstein 
distribution was imp leme nted in an analogous way as de- 
scribed in Appendix C 2 . 

As seen in Fig. [13, the transverse momentum depen- 
dence of the phase-space density is well described with 
our model whose geometrical and dynamical input was 
extracted from a space-time analysis of hadronic spec- 
tra. The chemical potentials needed to account for the 
absolute particle yields are listed for the three different 
parameter sets in Table |v|. In contrast to the box-shaped 
density profile, the multiplicities obtained from a Gaus- 
sian one require chemical potentials above 135 MeV. This 
may be taken as yet another hint that a box-shaped den- 
sity profile is favoured by the data. Also, at such high 
values of fj, the effects due to the finite fireball volume 
may become relevant and the treatment gets more com- 
plicated [Ell. 



3. Discussion 

The determination of the pionic chemical potential as 
a measure of the "overpopulation" of phase-space is of 
interest for the discussion of generic quantum- mechanical 
effects (e.g. stimulated emission [|5^]), for the study of 
the fireball's chemistry [ [53| , as well as for certain signals 
of in- medium changes of hadron masses [ ^4], |5^, Q . Here 
we discuss the implications of Fig. |l3| in this context. 

As noted above, the model curves in Fig. [l^ reproduce 
the transverse momentum dependence better than simple 
statistical distributions. However, they account for the 
absolute particle yields only with the help of a relatively 
large chemical potential while the statistical distributions 
rather overestimate the absolute yields without invoking 
a chemical potential. The latter effect can be understood 



1. The ^-dependence of the emission function (3.1) 
introduces the factor exp[— (77— 77o) 2 /(2(A?/) 2 )]. For 
the rapdity bin 3.9<y<4.4, this effectively amounts 
to a negative chemical potential of — 0.46T. For 
the temperatures from sets bl and b2 this gives 
w —50 MeV which must be compensated by the 
"usual" chemical potential. Note that this com- 
pensation would not be required if a box-shaped 
density profile were also used for the ^-distribution 
of the source; such a distribution seems also to be 
consistent with the recent study of the phase-space 
density in Pb+Pb collisions [ [l?J . 

2. The fitted chemical potential in Table |y| also ac- 
counts implicitly for the inclusion of pions from 
short-lived resonances. For high temperatures (~ 
150 MeV) this is an important effect, since then 
the contribution by these resonances is compara- 
ble with direct pion production. The exact fraction 
of particles from resonance decays depends on the 
chemical potentials for the resonances and requires 
a more detailed picture of the fireball chemistry. It 
decreases with decreasing temperature [ [57J . For a 
simple estimate one may count the number of pions 
from short-lived resonances as calculated from our 
model. The pion chemical potential needed to com- 
pensate the lack of these resonance contributions 
in our emission function by an additional amount 
of direct pions is fi « 20 MeV for the set bl and 
[i w 34 MeV for the set b2. The higher value of fi 
in the latter case reflects the larger resonance frac- 
tion at the higher temperature. 

Subtracting these two effects from the values given in 
Table we are led to a "real pion chemical potential" 
which is significantly smaller: we find n„ « 60 MeV for 
T = 100 MeV (set bl) and ^ m 25 MeV for T = 120 MeV 
(set b2). These parameters still indicate a slight overpop- 
ulation of the pion phase-space at freeze-out. Clearly, our 
discussion could be refined by a proper study of the ra- 
pidity dependence of the particle spectra and yields and 
by a more quantitative inclusion of contributions from 
short-lived resonances. This, however, lies outside the 
scope of the present work. 



VI. CONCLUSIONS 

In the context of simple analytical parametrizations 
of the freeze-out emission function, we have presented 
a space-time analysis of hadronic one- and two-particle 
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spectra measured for slightly forward rapidity by the 
NA49 Collaboration at the CERN SPS. Our work quotes 
for the first time full \ 2 confidence levels for the tempera- 
ture and the transverse flow. It allows for a quantitative 
discussion of possible model-dependencies via the com- 
parison of different transverse density profiles. Thus, our 
results corroborate quantitatively the picture of a colli- 
sion system with strong transverse collective flow which 
drives an expansion to twice the initial transverse size 
before freeze-out. During the evolution to this freeze-out 
stage the system seems to cool down significantly. 

Most importantly, this general picture does not depend 
on details of the model emission function with which 
the analysis is done. The \ 2 values of the fit are, how- 
ever, sensitive to details of the model. With the used 
transverse flow profile they favour a box-shaped over a 
Gaussian transverse density profile. This gives support 
to a similar conclusion reached recently in the study of 
deuteron coalescence models, and it indicates a possible 
different origin for negative (x 2 — y 2 ) contributions pre- 
viously attributed to opacity effects. 

The strong Mj_-dependence of the HBT side and out 
radius parameters requires a large transverse flow and 
thus favours a low kinetic freeze-out temperature in the 
range 80MeV < T < 110 Me V. Other earlier analy- 
ses which had obtained higher freeze-out temperatures 
therefore require some discussion: 

• In [~ 



] a temperature T — 120 MeV combined with 
v± = 0.43 was found from a simultaneous fit to 
single-particle mj_-spectra of different species. This 
value does not contradict our findings since it lies 
only slightly below the 90% confidence level of our 
fit to correlation radii. (The combination of trans- 
verse velocity and density profile used in [ |5^] dif- 
fers from ours. This can lead to slightly different 
fit results). While a simultaneous analysis of only 
the single-particle spectra of many different particle 
species thus also permits the separation of T and 
V_l, it must make the additional assumption that 
all these hadron species decouple simultaneously. 
That one obtains in this way similar results as by 
combining spectra and correlations of a single par- 
ticle species (pions) suggests that this assumption 
is in fact reasonably well justified. 

• By studying the curvature of the very accurately 
measured ir° spectrum in mj_ up to 4 GeV/c the 
WA98 collaboration extracted an apparently highly 
accurate value of the freeze-out temperature of 
about 185 MeV [|59|. This analysis was criticized [ 
|60] | on the basis that the curvature in the high p± 
region results from the Cronin effect and that an 
agreement with a thermal model is likely to be an 
artefact and is subject to severe systematic uncer- 
tainties in the model parametrization. 

• The NA49 collaboration has analyzed their data 
on the basis of simplified analytical expressions de- 



rived in [ [II]], such as e.g. Eq. (4.1). In this way 
they arrived at T » 120 MeV and v x » 0.55 [ [QJ. 
The discrepancy with our fit parameters can be 
traced back to the limited validity of these analyti- 
cal expressions, especially for the description of the 
h~ spectrum. If one compares their corresponding 
X 2 valley with ours (Figs. || and f|) one sees that 
it bends over at low T and large v± whereas ours 
doesn't. This is the main difference and accounts 
for their intersection to occur at somewhat larger 
values of T and v± . 



• An analysis in the spirit of ours was reported in [ 
Q. Again, higher temperatures (134-145 MeV de- 
pending on the analyzed set of data) were quoted. 
In this case we suspect that the high temperature 
values are driven by requiring a fit to the nor- 
malized single-particle spectra with vanishing pion 
chemical potential. The high temperature is thus 
a consequence of the observed multiplicity. More- 
over, the correlation radii were calculated from an- 
alytical approximations [ Q| which again introduce 
additional uncertainties. 



Our analysis of the observed phase-space density (5.16) 



has shown conclusively that the main deviation of the ex- 
perimental data from a purely quantum-statistical phase- 
space distribution is due to the collective dynamical ex- 
pansion of the system. Also, we have inferred from 
our analysis a pion chemical potential fi v ~ 60 MeV at 
T = 100 MeV. This indicates an overpopulation of the 
pion phase-space by about a factor 2 at thermal freeze- 
out, consistent [ |33[ with analyses of the particle ratios 
indicating a much earlier decoupling of the hadron yields 
(chemical freeze-out), namely at T c hcm ~ 170 MeV in 
Pb+Pb collisions at the SPS. Nevertheless, this pion ex- 
cess is too small to necessitate the study of stimulated 
pion emission [ B2| or higher order correlations [ k32| . 
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APPENDIX A: DISTRIBUTION OF THE PAIR 
MOMENTUM 

In this Appendix we derive distributions of the pair 
rapidity Y and of the pair transverse momentum K±. 
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The observed h single-particle rapidity spectrum can 
be very accurately approximated by a Gaussian distribu- 
tion of width Ay = 1.4 [ H, §7|: 



Qi(y) 



1 



exp 



{y - yo) 

' 2(A»)2 



(Al) 



Here yo stands for midrapidity, and the distribution is 
normalized to unity (rather than to the total number of 
particles) . The distribution of the pair rapidity 



Y = 



Vi + V2 



(A2) 



is then evaluated as 



1 



V* (AvY 



exp 



{Y-yof 
(Ay) 2 



(A3) 



When deriving the analogous distribution of the mean 
transverse momenta K± of the pairs one starts from the 
corresponding single-particle distribution (again normal- 
ized to unity) [ E7J 

here T; nv is the measured inverse slope, and Af is a nor- 
malization constant. The average transverse momentum 
of a pair with transverse momenta Pix,P2X is 

K± = \ (Pi± + Pa± + 2pi±p 2 ± cos(0i - 2 )) , (A5) 

where <f>\ and 02 are the azimuthal angles of the individ- 
ual transverse momenta. Then the desired distribution 
of K± is obtained from 

poo p2tv poo p2ir 

gi(Kx) = / dpix / #i / dp 2 ± / d(j} 2 
Jo Jo Jo Jo 

Xf>r(m) ei"(P2±) <5 (Aa - K x ) , (A6) 

AjL = \\jp\x+P\x + 2pi±P2± cos((/>i - 02 ) ■ (A7) 

We now integrate over the angles with the help of the 
(^-function. The integration boundaries are most easily 
implemented by using as integration variables 



0i_l+P2_l) =C e (K x , oo), (A8a) 
(pix - P2J.) = £ G {-Kx, Kx) • (A8b) 



Integrating over £ and substituting £ by a dimensionless 
variable a — C/Kx we obtain the final expression [ [l9| 

i K 3 , f°° , / 2if ± \ 2a 2 - 1 

g 2 {K ± ) = 2—± da exp -— - a . 

J i,w Jl \ J inv / Va - 1 



The resulting distribution, for an inverse slope T; nv = 
185 MeV as used in the calculations, is illustrated in 
Fig. It is compared with the two naive guesses (both 
normalized to unity) 



f{Kx) 
9{Kx) 



AKx 

w 

inv 

Kx 
T 2 

inv 



exp 



exp 



2K_ 
Ki 

1 irr 



(AlOa) 
(AlOb) 



Note that g{Kx) is actually the distribution of px inte- 
grated over the azimuthal angle (cf. ( A4)); it seems to ap- 
proximate Q^Kx) better than f(Kx)- However, as seen 
in Fig. 15 this is only the case for Kx < 500MeV/c. For 
high Kx the distribution g^Kx) asymptotically behaves 
like f(Kx)- Thus in general, a r eplac ement of g 2 {Kx) 
by either of the two distributions (A10) cannot be recom- 
mended. 



APPENDIX B: COMPUTATION OF THE 
NEGATIVE HADRON TRANSVERSE 
MOMENTUM SPECTRUM 

T he the oretical pj_-spectrum needed in the fit in Sec- 
tion 



IV A was calculated with the routine described in 



|37| . Here we point out the changes in the original routine 
which were made before it was used in this work. 



1. Bose-Einstein distribution for direct pions 

The directly produced pions are distributed according 
to the Bose-Einstein distribution, unlike in the original 
code in which the Boltzmann (high-energy) approxima- 
tion was employed. In practice this was done by ex- 
panding the Bose-Einstein distribution into powers of the 
Boltzmann distribution and truncating after a sufficient 
number of terms (see also Appendix |C 2| ). The pion chem- 
ical potential was set to zero. 



2. Chemical potentials for resonances 

Baryon number and strangeness conservation in the 
fireball requires the introduction of nonzero baryon and 
strangeness chemical potentials /ib and /is- Since these 
chemical potentials turn out to be not negligible f 



64 



, they may affect the multiplicities of bary- 
onic and/or strange resonances and thus also the number 
of pions. The condition of strangeness neutrality of the 
fireball determines /is as a function of T and [Ib- 

The yield of particles with strangeness S is then mul- 
tiplied by a factor As = exp(S//s(T)/T). Guided by the 
dependence of /is on the temperature shown in Fig. 3 of 
Rcf. [ p3| and with an eye on the results of [ p5|, H we 



(A9) introduced a polynomial parametrization for /is(T) 



12 



li B (T)^aT 2 + bT + c, 



(Bl) 



with a = 24.8 GeV , b = -7.94, and c = 0.905 GeV. Al- 
though rough, this parametrization is sufficiently precise 
for our purposes and easy to handle. The strangeness 
chemical potential fig was calculated from 

= ^5, ? ,™^e( fl -^+ s *)/ r K 2 (^) , (B2) 



where T is the freeze-out temperature and the sum runs 
over all species. Particle masses, baryon numbers and 
strangeness are denoted by rrii, Bi, and Si, respectively; 
<7i is the spin and isospin degeneracy. When determining 
Hs{T) all resonances up to 2 GeV were taken into ac- 
count, in order to be consistent with [ p3||. The resulting 
dependences of the chemical potentials on the tempera- 
ture are shown in Fig. O. 



3. Non-pion contributions and resonance decays 

Since we studied data on h~ -spectra rather than iden- 
tified pions, our calculations had to include non-pion neg- 
ative hadrons. Their distribution is assumed to be given 
by the same emission function with appropriately mod- 
ified masses and multiplied by the corresponding spin- 
isospin degeneracy factor. We used the Boltzmann ap- 
proximation, which is justified due to the large masses of 
these particles, with the appropriate chemical potentials. 

Some kaons and antiprotons are also produced by res- 
onance decays. We included the same set of resonances 
as in the calculation of the pion spectrum (see Table 1 
of [ p7|). The decays contributing to the K~ spectrum 
are li sted in Table [vj, those leading to antiprotons in 
Table [VT^ . 

As in [ p7y, contributions to the pion spectra from £ 
and A decays were treated as decays of a single resonance 
Y with an average mass of 1.15 GeV. Antiproton produc- 
tion via decay chains of the type Y,*^Y + ...^p+... 
was effectively included by enhanced branching ratios for 
the Y decay channels. (We have also replaced by J = 3/2 
the erroneous value J = 1/2 for the spin degeneracy of 
E* listed in [ [57).) As argued in [ [37), the above ap- 
proximations work well because they only concern small 
relative contributions to the full result. 



APPENDIX C: TOTAL PION MULTIPLICITY 

1. Boltzmann distribution 

In this Appendix we derive formulae for the total num- 
ber of produced pions predi cted by the model (3.1), 
both for the Gaussian (Eq. ( |3.3| )) and the box-shaped 



(Eq. (3.4)) transverse density profiles. Th ey are obtained 
by integrating the emission function ( |3.l| ) over positions 
and momenta: 



(N) = /o ^ 3 J pj^ dp±^ dy d(f> dip r dr di] 
(r-rof 



t dr 



(2tt 

x cosh(y — 77) exp 



2(At) 2 J 6XP V 2(At7) 2 
m± cosh(y - rj) cosh (77/ r/r rmB ) ' 



V24At) 2 
(V - Vo) 2 



m± 



x exp 



xG(r) exp 



T 

pj_ sinh(r) f r/r rms ) cos(<ft - <g) 
T 



The resulting expression takes the form [ [K 
<A) = JL TO 2 TK 2 (^) V inv , 



(CI) 



(C2) 



where Vi nv is the invariant volume of the fireball which 
takes into account the Lorentz-contraction of the moving 
volume elements and the Cooper-Frye-like flux through 
the freeze-out hypersurface [ Q . For a sharp freeze-out 
hypersurface 07 at r = const., for example, this leads to 



dV iri 



u • da , 



with 



d 3 a^ = (cosh rj, 0, 0, sinh 77) r dr\ r dr dip 



and (using (|3.5|)) 



dVw, 



u ■ d a = cosh rjtT dr\r dr dip . 



(C3) 



(C4) 



(C5) 



The model ( |3.1[ ) can be thought of as a Gaussian superpo- 
sition of r = const, hypersurfaces. Moreover, the volume 
is not populated equally densely and the corresponding 
distributions must also be taken into account: 



Wn 



2-7T j r dr cosh rj t (r) G(r) 

rdr ( (r-r ) 2 

exp 1 



^(At) 2 V 2 ( At ) 2 
(17- Vo) 2 ' 



dr] exp 



2(At7) 2 



(C6) 



The last two integrals are trivial, but different transverse 
density and flow profiles lead to different results for the 
first integral. The combination of fl3.6|) with (ETjl) yields 



' 1 



Gauss 
mv 



l( 27r ) 3r rmsn)A?7 



, i'Vf/2 



dx 



-rif/2 



(C7) 



while the box-shaped transverse density profile gives 



OX 

mv 



= (2^) 5 ^ ms r A?7 
v^sinh(v / 2r/ / ) cosh( v / 2r//) - 1 



(C8) 
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Here r rms = Rq\[2 for the Gaussian and r lma = Rb/V% 
for the box profile, respectively. For rjf — (no trans- 
verse flow) the box profile thus gives twice as many pions 
as the Gaussian one (at the same value of r rms , i.e. for 
the same interferometric signal). This is a volume effect: 
at the same r rms the invariant volume is twice as large 
for a box-shaped distribution than for a Gaussian one. 

Transverse flow increases the invari ant vol ume by the 
factor given in the square brackets of (C7 C8). The phy- 
sical picture is that the transversely moving fluid cells are 
Lorentz contracted and thus more of them are needed to 
"fill" the transverse profile of a given size. This Lorentz 
contraction is reflected by the factor coshry t (i.e. the 7- 
factor corresponding to the transverse motion) in (C5). 



2. Bose-Einstein distribution 

Since the pion mass is comparable to the freeze-out 
temperature, the use of the Boltzmann approximation 
for pions is questionable, in particular if they develop a 
positive chemical potential. In this subsection we show 
how the formulae for the total particle yield are modified 
if the Bose-Einstein distribution is taken into account. In 
this case the mean pion multiplicity is obtained from 



(It 



^2vr(A7 



exp 



\ / p ■ do 



exp 



p ■ u(x) 



T 



2(Ar) 2 
T 



(C9) 



The last integration gives the momentum spectrum emit- 
ted on a t = const, freeze-out hypersurface [ In the 
second integration one sums up contributions from such 
hypersurfaces with Gaussian distributed r values. The 
"usual" chemical potential is denoted by hq while the 
density distribution of the source is implemented via the 
"cc-dependent part of the chemical potential" [ ^3(| Q 
denoted by fJ.(x). The model with a Gaussian transverse 
density profile is thus characterized by 



p,{x) 



(v - m) 2 



2R% 



T 2(At7) 2 
while the box-shaped source is implemented via 



jti(x) (77 - r/ ) 2 



2(At7) 2 



B(r), 



where 



B(r) 



if r < R B , 
00 otherwise. 



(CIO) 



(Cll) 



(C12) 



For practical evaluation the Bose-Einstein term under 
the last integral in (C9) is expanded into a geometric 



series. When exchanging the order of summation and 
integrations one arrives at 



d 3 p 



dr 



V2tt A? 



exp 



p ■ d a exp 



p ■ u(x) - no 
(T/n) 



2(At) 2 
(T/n) 



(C13) 



The Boltzmann approximation used in the previous sub- 
section is recovered from the first term n = 1. Higher 
terms n > 1 can be regarded as Boltzmann contributions 
with lower effective temperatures ~. In practice, for suf- 
ficiently small chemical potentials, this allows to truncate 
the expansion at a certain value nt rU nc since the size of 
the contributions decreases exponentially with n. For a 
Gaussian density distribution one finds explicitly 



(TV) 



T o r 2 ms A77 m 2 T n 

n=l 



5/2 e »»/T K « 



Vf 

2Vn 



dx 



\ T J 



(C14) 



Note that the invariant volume cannot be factorized from 
the expression. This is generally true for models with a 
Bose-Einstein distribution with an x-dependent chemical 
potential. A special case is the box-model where the 
invariant volume (C8) still factorizes: 



(AO = V, 



box 



m 2 T 



71=1 



yj, -3 / 2 



/r K 2 f^). (C15) 



\ T J 



The different factorization properties of (C14) and ( C15| ) 
imply that for 77/ = the yields are no longer related by a 
simple factor 2 once Bose statis tics is taken into account. 

Expressions ( |C14 ) and (C15) were used in_the calcula- 
tions of the direct pion multiplicity in Sec. [V C . There, 
the chemical potential (io was set to since che mical 
equilibrium was assumed. The calculations in Sec. VB 2 
were done with the values of fj, = fj,o given there. 
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TABLE I. Correlation radii for the Cartesian parametriza- 
tion in the LCMS obtained by calculating the averages of 
the squared radii from the analysis of MTPC and VTPC 
data (see text). In the second row the values of K^ n for 
each bin according to ( [2.2; ) are displayed. Selected rapid- 
ity bin: 3.9 < Y < 4.4. (In fact, in the MTPC analysis 
the rapidity bin 4 < Y < 4.5 was used. Here that from 
VTPC data sets has been adopted.) Momenta are given 
in units of GeV/c. In the K_i_-bin 0.3 — 0.45, data from 
[ H for 0.3GeV/c < K ± < 0.4GeV/c and from [ || for 
0.3 GeV/c < K± < 0.45 GeV/c were averaged. 



K±-bin 


- 0.1 


0.1 - 0.2 


0.2 - 0.3 


0.3 - 0.45 


K bin 


0.065 


0.151 


0.247 


0.365 


Rs (fm 2 ) 


33.3 ±4.0 


27.0 ±3.5 


25.0 ± 1.8 


21.4 ±2.5 


Rl (fm 2 ) 


38.7 ±2.1 


38.8 ± 3.3 


33.6 ±4.5 


30.4 ±5.2 


Rf (fm 2 ) 


61.9 ±5.3 


46.7 ±5.0 


27.4 ±5.6 


20.7 ±3.1 


Rli (fm 2 ) 


8.4 ±3.0 


13.8 ± 2.2 


8.4 ±3.7 


7.2 ±5.6 



TABLE II. The same as Table 02, but for the correla- 



tion parameters in the Yano-Koonin-Podgoretskii parametri- 
zation. v YK is measured relative to the LCMS. The same 



comments about binning as in Table C 2 apply 



K±-bm 


- 0.1 


0.1 - 0.2 


0.2 - 0.3 


0.3 - 0.45 


K bin 


0.065 


0.151 


0.247 


0.365 


Rs (fm 2 ) 


32.7 ±4.7 


28.4 ±3.3 


26.0 ±3.6 


23.2 ±5.5 


Rl (fm 2 ) 


59.6 ±9.8 


42.3 ± 5.8 


25.9 ±6.7 


21.1 ±3.3 


Rl (fm 2 ) 


28.6 ± 18.7 


14.8 ± 9.2 


8.4 ±6.7 


8.1 ±5.7 




-.15 ± .04 


-.26 ± .04 


-.23 ± .05 


-.25 ± .10 



TABLE III. Sets of model parameters for box-shaped (b*) 
and Gaussian (g*) transverse density profiles, as obtained 
from fits to the single-particle mx-spectrum and two-particle 
correlation radii. For each parameter set, the average trans- 
verse expansion velocity v± and the total n~ multiplicity N 
(cf. Fig. 0) was calculated from the model. The errors for 
N include only the uncertainties in the values of the model 
parameters. 



Box set 


bl 


b2 


b3 


T (MeV) 


100 


120 


160 


Vf 


0.6 


0.5 


0.35 


Rb (fm) 


12.1 ±0.2 


11.5 ±0.2 


10.7 ±0.2 


To (fm/c) 


6.3 ± 1.1 


5.5 ±1.1 


4.4 ±3.5 


At (fm/c) 


3.6 ±0.6 


3.2 ±0.7 


2.6 ±2.0 


Ar/ (fixed) 


1.3 


1.3 


1.3 


Vj_ 


0.5 


0.43 


0.33 


N 


139 ± 24 


254 ± 56 


813 ± 649 


Gauss set 


gl 


g2 


g3 


T (MeV) 


100 


120 


160 


Vf 


0.6 


0.48 


0.35 


Rg (fm) 


6.5 ±0.1 


5.9 ±0.2 


5.6 ±0.1 


t (fm/c) 


7.8 ±0.8 


6.6 ±0.9 


5.5 ±0.9 


Ar (fm/c) 


2.3 ±0.7 


2.3 ±0.7 


1.8 ±0.8 


Ar/ (fixed) 


1.3 


1.3 


1.3 


Vj_ 


0.46 


0.39 


0.29 


N 


96 ± 10 


157 ± 24 


544 ±91 



TABLE IV. The values of the effective intercept parameter 
as a function of K±, estimated from the NA49 data. 



K± (GeV/c) 


0.065 


0.151 


0.247 


0.365 




0.4 


0.495 


0.46 


0.4 



TABLE V. Estimated chemical potentials needed for the 
box-models of Table |ffl| to fit the data for the average 
phase-space density. 



set 


T (MeV) 


fj, (MeV) 


bl 


100 


123 ±9 


b2 


120 


115 ± 15 


b3 


160 


90 ±20 



1G 



TABLE VI. Decay channels included in the calculation 
of the K~ contribution to the h~ -spectrum. Branching ra- 
tios were multiplied by the Clebsch-Gordan coefficients cor- 
responding to the given isospin states. 



channel M res0 n (MeV) 


T (MeV) 


J 


BR 


K*~ -» n u K~ 892 


115 


1 


1/3 x 1 


K*° -» tt+K- 892 


115 


1 


2/3 x 1 


TABLE VII. Decays contributing to the product 


ion of an- 


tiprotons. Branching ratios were multiplied by the appropri- 


ate Clebsch-Gordan coefficients. 








channel M Icson (MeV) 


r (MeV) 


J 


BR 


A++ _> Tr-p 1232 


115 


3/2 


1 


A+ -> ty°P 1232 


115 


3/2 


2/3 x 1 


A -> ty+p 1232 


115 


3/2 


1/3 x 1 


A^vr+p 1116 


» 


1/2 


0.639 


E+ — ► ty°p 1193 


» 


1/2 


0.516 


E° -> A7 -> 7T+p 1193 


» 


1/2 


0.639 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 

K ± (GeV/c) K ± (GeV/c) 

FIG. 1. The correlation radii in the Cartesian parametriza- 
tion in the LCMS (see text for d etai ls) . Black symbols indi- 
cate the data points from Table |C 2j . Gray dashed symbols 
represent the values calculated with cross-check rel ation s from 
the measured YKP parameters given in the Table C 2 . 



04 o 



E 

CM = 



0.1 

K, 



0.2 0.3 
(GeV/c) 



0.4 



FIG. 2. The measured YKP parameters 




given 



in the LCMS. Black symbols indicate the data points from 
Gray dashed symbols represent the values calcu- 



Table C2 



lated with cross-check relations from the measured Cartesian 



parameters given in Table C 2 
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0.5 






0.4 Vl 
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0.3 






0.3 


0.2 






0.2 



0.08 0.1 0.12 0.14 0.16 
T (GeV) 



(b) 

FIG. 3. x 2 contour plots for the fit of the single-particle 
h~ spectrum to the model ([Tl]) with a box-shaped transverse 
density profile. The contour lines correspond to the following 
confidence levels: 99%, 90%, 50%, 10%, 5%, 1%, 0.5%, 0.1%, 
0.05%, and 0.01%. Panel (a): fit including all pions. Panel 
(b): negative hadrons, including also K~ and p and their 
resonance decay contributions. 



1F 




2- GAUSS, all pions 

I I I I I L_ 

0.08 0.1 0.12 0.14 0.16 
T (GeV) 



(a) 

1F 




2-GAUSS, all h" 

i i i i i i 

0.08 0.1 0.12 0.14 0.16 
T (GeV) 

(b) 

FIG. 4. The \ 2 contour plots as in Fig. fcflbut for the model 
with Gaussian transverse density profile. The average 
transverse expansion velocity v± allows for a direct compar- 
ison between the two models. The corresponding values of 
rjf are also shown; they differ substantially between the two 
models. 
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FIG. 5. x~ 

to the model ( p.f| ) with a box-shaped transverse density pro- 
file (solid lines). The numbers give the corresponding confi- 
dence levels. The dashed lines show the result of the fit to the 
single-particle spectrum (Fig. |^b). Crosses denote the posi- 
tions of the three parameter sets used in further calculations. 




0.3 



0.2 



^10% v \VoV;\\ \NX 



GAUSS 



0.4 v ± 



0.2 



0.08 0.1 0.12 0.14 0.16 
T (GeV) 

FIG. 6. Same as Fig. B but for the model with a Gaussian 
transverse density profile. 




0.4 Vl 



0.08 0.1 0.12 0.14 0.1 
T (GeV) 



0.08 0.1 0.12 0.14 0.16 
T (GeV) 



FIG. 7. Total ir~ multiplicity (resonance decays included) 
as a function of T and v±, computed from the model with 
(a) box-shaped, (b) Gaussian transverse density profile. The 
wiggles in the contours towards the right in panel (a) result 
from numerical uncertainties in the fitted model parameters. 
For other details see text. 



FIG. 8. Measured h spectrum in the rapidity window 
4.15<t/<4.65 [ ^ij compared to the best model fits bl and 
gl- 




0.2 0.3 0.4 
M x [GeV/c 2 ] 

FIG. 9. Comparison of model predictions from the model 
with box-shaped transverse geometric profile and inserted pa- 
rameter sets bl, b2, b3 with the data on Cartesian correlation 
radii (to which the model was fitted). 




FIG. 10. Same as Fig. ^ but for the Gaussian transverse 
density profile with parameter sets gl, g2, g3. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

m ± - m„ [GeV/c 2 ] 




o 

x (fm) x (fm) 

FIG. 11. Transverse cuts through the effective emission re- 
gion of particles with _Kx=500 MeV/c. Left: Gaussian trans- 
verse density profile. Right: box-like transverse density pro- 
file. In both cases T=100MeV, ??/=0.8. 
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FIG. 12. Schematic sketch of the interferometric lifetime 
measurement. The time parameter to is extracted under the 
assumption that the longitudinal expansion is boost-invariant 
during the entire dynamical evolution. If the boost-invariant 
velocity profile develops gradually, the fireball can live longer. 
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FIG. 14. The distribution g 2 (K ±) of average transverse 
pair momenta calculated from (|A9|). The dashed line shows 
the function f( K±) fr om Eq. ( |A10a| ) ; the dotted line stands 
for g(K±) from ( |A10b| ). T inv = 185 MeV. Longer tics indicate 
the bin positions for the data in Sec. [lj|. 




P ± [GeV/c] 



0.1 0.2 0.3 0.4 

P ± [GeV/c] 

FIG. 13. Average phase-space density (f){p) in the rapid- 
ity bin 3.9 < Y < 4.4 as a function of p±. The data points are 
determined from NA49 data (Sec. VB 1). Black lines d enote 
expe ctatio ns from the statistic al tre atment in Sec. V B 2: (/s) 
(eajlH), solid); (f B ,) (eq. ( p^ ), long-dashed); (f BE ) (eq. 
(5.18), short-dashed); and {/be'} (eq- (5-20), dash-dotted). 



Results of our model (see text) are shown by thick solid gray 
lines. The symbols in the upper right corners indicate the used 
set of model parameters; the corresponding pionic chemical 
potentials are listed in Table [v| The right column gives a 
logarithmic representation of the left column. 




0.6 0.8 1 
K, [GeV/c] 

FIG. 15. The same as Fig. [14l but all three distributions 
are divided by K± and plotted logarithmically. 




0.06 0.08 0.1 0.12 0.14 0.16 0.18 
T [GeV] 

FIG. 16. The chemical potentials (j,b(T) and fJ,s(T) in units 
of T as given by Eqs. (M) and (Bgh. 
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